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ABSTRACT 

• In the standard ACDM cosmological model with a Gaussian primordial density fluc- 

""q^I tuation field, the relatively low value of the mass variance parameter (o-g = 0.74lQ Qg, 

I ■ obtained from the WMAP 3-year data) results in a reduced likelihood that the mea- 

O I sured level of CMB anisotropy on the scales of clusters is due to the Sunyaev-Zeldovich 

^ (S-Z) effect. To assess the feasibility of producing higher levels of S-Z power, we explore 

I two alternative models which predict higher cluster abundance. In the first model the 

primordial density field has a Xi distribution, whereas in the second an early dark 
energy component gives rise to the desired higher cluster abundance. We carry out 
^ i the necessary detailed calculations of the levels of S-Z power spectra, cluster number 

' counts, and angular 2-point correlation function of clusters, and compare (in a self- 

. consistent way) their predicted redshift distributions. Our results provide a sufficient 

I basis upon which the viability of the three models may be tested by future high quality 

measurements. 

' Key words: galaxies: clusters: general - cosmic microwave background - large-scale 

, structure of the Universe 
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^ 1 INTRODUCTION 

C3 Measurements with the BIMA (Dawson et al. 2002), CBI (Readhead at al. 2004), and ACBAR (Kuo et al. 2004) experiments 
indicate a significant power excess at high multipoles (£ > 5000 in the former, I ~ 3000 in the two latter cases) with respect 
to theoretically predicted levels of the primary CMB anisotropy. It has been suggested that the excess could be due to the 
S-Z effect, but its observed level would require a rather high value of ag (> 1) in order to have the necessary large number of 
massive clusters, with the largest relative contribution to the S-Z power spectrum. Such a high value is usually disfavoured 
in cluster studies, and is at clear variance with results from the latest WMAP 3-year data, according to which erg = 0.74^Q Qg 
(Spergel et al. 2006). If the reported excess is indeed due to the S-Z effect, modifications to the standard ACDM cosmological 
model may be required in order to account for the implied higher number of clusters. 

Cluster abundance at high redshift can be boosted by the presence of a scale-dependent, non-Gaussian, positively skewed 
component in the primordial density fluctuation field, as was demonstrated by Mathis, Diego, & Silk (2004), and Sadeh, 
Rephaeli, & Silk (2006) (hereafter SRS). In such a model, primordial overdensities attain larger amplitudes with higher 
probabilities than the corresponding overdensities in a purely Gaussian random field, and may therefore give rise to earlier 
cluster formation and to higher numbers of massive clusters. In both works the former WMAP 1-year normalisation was 
employed {og, = 0.9), and it was shown that already with this relatively high value it is difficult to reconcile the CMB power 
excess with the inferred cluster population, if the underlying fiuctuation field obeys Gaussian statistics. On the other hand, 
a Xi distributed field was demonstrated to be capable of producing a sufficiently large cluster population, and consequently 
S-Z power levels that are consistent with the above mentioned observational results. The Xm family of models was originally 
proposed by Peebles (1997, 1999a, b), and is based on an isocurvature cold dark matter (CDM) scenario in which the primordial 
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density field is proportioinal to the square of a random Gaussian process (Peebles 1997). Implications of this model on the 
primary CMB anisotropy and the dynamics and evolution of the large scale strcuture have been the subject of several studies 
(e.g. Koyama, Soda & Taruya 1999, Mathis, Diego & Silk 2004). Other theoretically motivated non-Gaussian models found in 
the literature include, e.g., the log-normal probability distribution function. A rather comprehensive list of such models and 
their statistical properties are explored in detail in Coles & Barrow (1987). A common parametrization of non-Gaussianity is 
provided by the following transformation on an underlying Gaussian field: 

^{x)=acl>{x) + e[ct>^{x)-{ct>^}]. (1) 

The parameter e, alternatively written as /nl and referred to as the nonlinear coupling factor, characterises the amplitude of 

the quadratic term. In this respect, the Xm model constitucs a special case in which a = 0. The degree of non-Gaussianity of a 
given model is defined in several ways; these include the skewness of the density field (for which a one-to-one correspondence 
with the nonlinear coupling term can be derived; Matarrese, Verde & Jimenez 2000), and the parameter ^ defined as 

which describes the 3-a tail excess probability of the non-Gaussian model over the Gaussian probability distribution function. 

For the Xi model C, = 16.3. While there seems to be observational evidence for a lower level of non-Gaussianity (c^ < 4, 
Robinson & Baker 2000; Avila-Reese et al. 2003, and C ^ 1-6, Robinson, Gawiser, & Silk 2000; Komatsu et al. 2003), we 
believe that these results do not necessarily rule out models with higher levels of non-Gaussianity. This is due to the fact that 
these results were inferred either from primary CMB anisotropy observations, which explore scales much larger than those 
associated with clusters of galajcies, and therefore, lower overdensities, which may be indistinguishable from a non-Gaussian 
distribution in case the density field is scale- variant, or from observations of low-redshift clusters, whereas - as will be shown 
below - the non-Gaussian tail would be mainly manifested in an enhanced high-redshift cluster population. 

Early dark energy (hereafter EDE) models provide an alternative for generating higher cluster populations at higher 
redshifts. In these models the dark energy density is appreciable already at early epochs of cosmic evolution, and attains the 
observationally inferred value at present. Consequently, the quintessence equation of state coefficient changes with time. Such 
models have been the subject of study by several authors who investigated their potential infiuence on the CMB (Doran et al. 
2001a, Caldwell et al. 2003), and the linear growth of structure (Ferreira & Joyce 1998, Doran et al. 2001b). More recently, 
Bartelmann, Doran, & Wetterich (2006; hereafter BDW) have carried out a comprehensive study of two specific EDE models, 
evaluating numerically the quantities relevant to structure formation and the halo mass function, such as the linear growth 
factor of density perturbations, critical density for spherical collapse, and the overdensity at virialisation. Specifically, they 
find that with respect to their corresponding values in the standard ACDM model, a slower evolution of the linear growth 
factor and reduced values of the critical density for spherical collapse are predicted. As noted by BDW, the slower evolution 
of the linear growth factor in EDE models is a consequence of the higher expansion rate of the universe at early times due to 
an early acceleration phase caused by the non-vanishing dark energy component. Therefore, with a given value of the present 
mass variance normalisation the corresponding quantity at early times should be larger than what is implied in the ACDM 
model. Since the critical overdensity for collapse at a given redshift is the linearly extrapolated value from the early universe, 
the slower evolution of the growth factor in EDE models is manifested by a lower critical density with respect to the ACDM 
model at all redshifts relevant to structure formation. Hence, the cluster population, which increases with decreasing 5c and 
increasing ctm, grows considerably. 

In this paper we focus on the impact of these two alternative scenarios on S-Z observables, namely the angular power 
spectrum of CMB temperature variations due to the S-Z effect, and the angular 2-point correlation function (hereafter A2PCF) 
of clusters. As is demonstrated below, the higher abundance of massive clusters at high redshifts in the non-Gaussian and 
EDE models boosts S-Z power levels (as a result of higher values of the Comptonization parameter), which in turn gives rise 
to larger temperature variations in the CMB. Conversely, the A2PCF is demonstrated to predict lower correlation levels in 
the EDE, and in particular, the non-Gaussian model, owing to the detailed properties of the linear bias parameter, which, 
with respect to the ACDM model, has a different mass and redshift dependence in the EDE model, and a distinct functional 
form in the nori- Gaussian model. 

In §2 we briefly describe the Press & Schechter mass function variant adapted to the Xi-distributed probability distribution 
function (PDF) of the primordial density fluctuation field, and outline the properties of the EDE model adopted in our 
calculations. Results of the calculations are described in §3, where S-Z power levels calculated for the standard ACDM, EDE, 
and non-Gaussian models are compared with the BIMA, CBI, and ACBAR observational results, and an account of the 
A2PCF for these 3 models is given. §4 includes a discussion and our conclusions. In the appendix wc provide a brief exposition 
of the differential equations governing the redshift evolution of the comoving radial distance, linear growth factor of density 
perturbations, critical density for spherical collapse, and overdensity at virialisation corresponding to the EDE model, and 
present some numerical results. 
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n(M, z) = ^^,F{^)^y^^—dM, (3) 



2 METHOD 

The standard ACDM, EDE, and non-Gaussian models will be referred to hereafter as models I, II, and III, respectively. A 
full account of models I and III is provided in SRS; here we only include new relevant aspects of these two models. We adopt 
a Press & Schechter (1974) mass function of the form 

pb da 
' Malj(z) dM 

where /i = 5c{z) /aM^z) is the critical overdensity for collapse in terms of the mass variance ctm at redshift z, and pb is the 
background density at z = 0. For models I and II, 

F{p) = \[^e~^-'''\ (4) 

whereas the corresponding expression for model III, for which we take a primordial density fluctuation field obeying Xi 
statistics, is 

(l+y2M)/2 

F{p) = , (5) 

erfc[l/\/2]y7r(l + V2Ai) 

The factor erfc[l/v'2] (with erfc representing the complementary error function) in the denominator is introduced in order to 
arrange that all the mass in the universe be included in haloes, as is similarly accomplished by introducing the well known 
factor of 2 in the original Press & Schechter mass function. (Recall that the full theoretical justification for this factor is 
somewhat uncertain; here it is merely introduced to ensure consistency with the mass function of model I). The predicted 
statistical properties of S-Z observables (i.e. weighted by the cluster population) will be affected by the model of the mass 
function; clearly, model III differs from models I and II in the functional form of the mass function. Model II differs from 
both models I and III in the redshift dependence of the critical density for collapse and the linear growth factor of density 
perturbations, as will be demonstrated below. 

For the calculation of the mass variance or we used a top-hat window function and CDM transfer functions taken from 
Bardeen et al. (1986): adiabatic transfer function for models I and II, and isocurvature transfer function for model III. Note 
that the presence of EDE is likely to affect the shape of the CDM transfer function; this effect is minor, however, with a 
largest difference of < 7% at the lowest wave numbers (fc ~ 10~*Mpc~^/i), becoming negligible at higher wave numbers 
(which correspond to cluster scales), as was verified using the latest version of the CMBFAST code (Zaldarriaga & Seljak 
2000), which admits a time-dependent equation of state for dark energy. Consequently, it is reasonable to use the fit provided 
by Bardeen et al. even in the case of model II. 

The calculations were carried out using cosmological parameters deduced from the WMAP 3-year data: Ha = 0.76 (which, 
for the case ui / — 1 is usually written as 0.q), 57™, = 0.24, h — 0.73, and as = 0.741q;o6- In addition to these parameters we 
take the spectral index to be n = 1 in models I and II, and n = —1.8 in model III. Model II requires additional parametrization 
of the EDE. We use a variant of the EDE models studied in BDW, where we take the density of early quintessence with 
respect to the critical density to be Qe = .0008, and the equation of state coefficient at z = to be mo = —0.99. The effective 
coefficient as a function of redshift is given in the expression (Wetterich 2004) 

= i+.io7(i+.) ' 
"^^'-^"" .o.(i^y:Z(i^) - 

The differential equations governing the evolution of the linear growth factor, critical density for spherical collapse, and 
overdensity at virialisation pertaining to this model, as well as their numerical solutions, are presented in the Appendix. It 
is important to note that while our choice of parameters in model I is completely consistent with WMAP results, this is not 
necessarily the case in models II and III. Nevertheless, we believe that we are able to account for this apparent inconsistency 
by choosing to normalise the resulting mass functions so as to yield the same (cumulative) cluster density at z = 0, as 
predicted by model I. This is justified given the fact that the observables quantified in this study are all related to clusters 
of galaxies. In fact, it turns out that the same normalisation of as = 0.74 (in accordance with WMAP results) induces an 
identical cumulative cluster density population at z = 0, as can be clearly seen by inspection Fig. [1] We include for reference 
a fourth model, which is basically identical to model I, but with a different normalisation of as = 0.8. This higher value 
has been advocated by Evrard et al. (2007), who claim that with a low normalisation, the specific energy in the IC gas and 
galaxies must increase by 50%, which would imply high biases with respect to the cluster dark matter, at variance with what 
is observed in hydrodynamical simulations. According to Evrard et al., a higher normalisation can solve this conflict, and also 
provides a better match to recent S-Z observations. These results are actually in line with several cluster studies which seem 
to indicate higher normalisations than inferred from CMB observations. 
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The S-Z angular power spectrum was calculated using a formalism similar to that adopted by (e.g.) Cooray (2000), 
Komatsu & Seljak (2002), and Refregier & Teyssier (2002). Note that these authors have used the NFW formalism to 
characterise the cluster dark matter profile and the implied IC gas profile. In our calculations we chose to model the IC gas 
with an isothermal /3-King profile with temperatures derived within the framework of hydrostatic equilibrium. The reason 
for our selection stems from the fact that the redshift and mass modeling of the concentration parameter appearing in the 
NFW profile have been carried out on the basis of results from N-body simulations, which usually incorporate primordial 
Gaussian fields. In such simulations clusters show up only at relatively low redshifts, such that it is not clear whether the 
deduced fitting of the concentration parameter can really address also higher-redshift clusters, a natural outcome of positively 
skewed non-Gaussian models. Furthermore, fitting formulae of c(M, z) found in the literature fail to yield meaningful values 
at redshifts at which clusters may already form in non-Gaussian models. Calculation of the A2PCF of (S-Z) clusters follows 
the work of SRS and that of Mei & Bartlett (2003); full details can be found in these papers. 



3 RESULTS 

Before we present results of the S-Z power spectrum and A2PCF, it is useful to investigate the behaviour of the mass function 
in the 3 models at various redshifts. Fig. [T] depicts the cumulative mass function of clusters with masses lying in the range 
lO"M0/i-^ ^ M s: IO^^Mq/i"^ at redshifts z = 0.01 and z ^ 3. As can be clearly seen, the three models are correctly 
normalised so as to yield the same cumulative cluster density at low redshifts. As expected, models II and III consistently 
predict higher abundances of high-mass clusters than the corresponding high-mass population in the standard ACDM model, 
with model III dominating over model II. This occurs at both low and high redshifts, and becomes more pronounced with 
increasing redshift. It appears then that the presence of an excess of high overdensity fluctuations associated with the Xi 
model has a stronger impact on the population of massive clusters than the slower evolution of the linear growth factor and 
lower values of the critical density for collapse that characterise model II. 

The behaviour of the respective mass functions is directly reflected in the S-Z power spectra illustrated in Fig. [2l Here 
the enhanced massive cluster population, particularly at relatively high redshifts, is manifested in both increased power levels 
(~ 1 • 10~^^, 5 • 10~^^, 3 • 10"^" in models I, II, and III, respectively), and a shift of the peak power towards higher multipoles 
from £ ~ 4000 for models I and II, to ~ 7000 for model III, reflecting the higher abundance of distant, low angular size 
clusters, particularly in model III. It can also be seen that model IV predicts power levels which are a factor 1.8 —- 2 higher 
than those obtained for model I. The numerical results for the four models were scaled to a frequency of 31 GHz (at which 
the non-relativistic spectral distortion due to the thermal S-Z effect amounts to a multiplicative spectral factor of ~ 3.8 in 
the expression for the power spectrum) so as to correspond to BIMA, CBI, and ACBAR observational results (that are also) 
shown in the figure. Of the four models, model III provides the best match to the observed CMB power excess, particularly so 
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Figure 2. S-Z angular power spectrum at u = 31 GHz in models I (continuous), II (dashed), III (dash-dotted), and IV (thick dash- 
dotted). The shaded areas correspond to the WMAP reported l-cr errors in eg. Also shown is the power excess reported by the BIMA 
(diamonds), CBI (crosses) and ACBAR (x-symbols) experiments. Note that the BIMA indicated power at £ = 8748 is an upper limit. 



at £ ~ 1000 — 1500 (at higher multipoles the calculated power levels are actually higher than the CBI and ACBAR results), 
whereas the predicted power in model 1 is far below and seems clearly inconsistent with the observed level. Although the 
predicted power levels in model II are higher than in model I, these are still only marginally consistent with the observational 
results. 

While both the EDE and non-Gaussian models generate higher S-Z power levels with respect to those predicted by the 
standard ACDM model, the levels of the A2PCF of S-Z clusters are lower in model II, and particularly so in model 111, with 
respect to the corresponding levels predicted by model 1, as illustrated in Fig. [3] The A2PCF of models 1, 11, and III peak 
at w ~ 2, 1, 0.5 at an angular separation of 1', falling off towards levels of w ~ 0.4, 0.2, 0.001 at angular separation 100'. Also 
noticeable is the fact that the curve associated with model 111 has a greater slope than the corresponding slopes of models I 
and 11. In order to understand these properties, it is useful to consider the redshift and mass dependence of the linear bias 
factor, which in the Gaussian case assumes the form (Mo & White 1996): 

KM,.) = l + t^, (7) 

where ^ = Sc{z)/[aM{0)D+{z)], D+{z) denotes the linear growth factor of density perturbations, and aA/(0) is the mass vari- 
ance at z = 0. Koyama, Soda, & Taruya (1999) generalised this expression to also address non-Gaussian models. Implementing 
their expression to the Xi model, we have 

6(M, z) = l + . (8) 

The lower bias levels of massive haloes in non- Gaussian models are a consequence of the linear bias theory developed by Mo 
& White (1996), and were explored in detail by Amara & Refregier (2004) within the framework of the log-normal PDF. In 
fact, these authors showed that for a log-normal PDF with an excess probability in the 3ct tail of ~ 10, the corresponding 
bias factor is larger than the one predicted by the Gaussian model at halo masses > 5 • 10^^ Mq h-^. Since the non-Gaussian 
model considered here has ~ 16, the mass scale at which the bias factor begins to be dominated by the corresponding 
Gaussian bias reduces to ~ 10^"^ AIq , which is the lower limit of our mass integral. It is therefore safe to assume that the 
relevant bias factor is lower than that of the Gaussian model within the entire mass range. 

The properties of the A2PCF can be now easily explained by inspection of the functional dependence of the bias factor 
on mass and redshift. The bias factor increases with increasing mass, since the mass variance decreases monotonically with 
increasing mass. Additionally, the bias factor increases with increasing redshift since the critical density for collapse increases 
with increasing redshift, whereas the linear growth factor decreases with increasing redshift. These properties are common 
to all of the three models. The differences between the predictions of the 3 models arise either from the different evolution 
of Sc{z) and D{z) in model II, or the distinct functional form of the bias factor in model 111. In model II, the magnitude of 
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Figure 3. The A2PCF of clusters calculated from their S-Z fluxes in models I (continuous), II (dashed), and III (dash-dotted), as a 
function of separation angle, 9. The calculations were carried out for the 353 GHz channel of the Planck/HFI experiment, with a beamsize 
of 5.1' and a limiting flux of 30mJy. 

the bias factor is lower than in model I, in hght of the lower and higher values of the critical density for collapse and linear 
growth factor, respectively (Fig. [5]). This is reflected in the lower levels of the A2PCF in model II, with respect to model I. 
More striking are the differences between the predicitons of models I and III. As can be easily realized from equation ((8]l , the 
bias is reduced owing to the addition of a larger than unity factor in the denominator of the second term, which is reflected 
in the substantially lower levels of the A2PCF. Note however, that the differences between the correlation levels of models I 
and III increase with increasing angular separation. At low angular separations both nearby and distant clusters contribute to 
the A2PCF by virtue of their relatively low mutual distances. Since the occurence of high-redshift, massive clusters, capable 
of generating the required flux limit, is substantially larger in model III, these contribute signiflcanly to the correlation levels 
of model III at low angular separations, whereas the corresponding contribution in model I originates in clusters residing 
at lower redshifts. At higher angular separations, where correlations can only arise locally, a relative deficit of low redshift 
clusters appear in model III, as can be seen in Fig. |4l where we plot cumulative number counts of S-Z clusters (evaluated at 
a frequency of 353 GHz, with beam size of 5.1', and flux detection limit 30 mjy). Indeed, if we compare the upper left-hand 
panel with the lower panel of Fig. 14] we may see that while (cumulative) cluster counts continue to increase in the redshift 
range 0.2 ^ z > 0.01 in model I, they level off at z < 0.2 in model III. This property of model III gives rise to the faster 
descent of the A2PCF, with respect to model I, at high angular separations. 

The results presented so far focused on the differences in S-Z power and angular correlation levels in the three models 
due to the respective mass functions and bias factors. Internal cluster properties relevant to the S-Z effect, such as virial size 
and the temperature and density of the intracluster gas, may be affected as well, as can be seen in the bottom panel of Fig. O 
where A^, (z), the overdensity (scaled to the critical density) at virialisation is plotted as function of redshift. However, the 
variation of A„(z) among the models amounts to less than 4%, implying little change in the above cluster internal properties, 
particularly so in comparison with the major differences apparent in the mass function and bias parameter. 



4 DISCUSSION 

The reported low value of as inferred from the WMAP 3-year data implies a signiflcantly reduced population of high-mass 
clusters of galaxies, as predicted by a Press & Schechter mass function based on a Gaussian density fluctuation fleld. With this 
low value it becomes even more difficult to attribute the CMB power excess measured with the CBI and ACBAR experiments 
at high multipoles to the (thermal) S-Z effect. The reason for this difficulty lies in the fact that the largest contribution to 
S-Z power comes from high-mass clusters, whose number density decreases steeply with decreasing as- If indeed this excess 
is attributed to the S-Z effect, an enhanced massive cluster population at moderate and high redshifts is required. In this 
work we have studied the implications of two non-standard cosmological models on S-Z observables, specifically, the angular 
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Figure 4. Cumulative number counts of S-Z clusters with fluxes exceeding 30 mjy as a function of redshift. Tiie calculations were 
carried out at a frequency of 353 GHz and beam size 5.1'. Upper left- and right-hand panels and the bottom panel correspond to models 
I, II, and III, respectively. Continuous curves depict the contribution to the cumulative number counts in the entire mass range of 
lO^^M0h~^ ^ Af ^ 10^^ MQh~^ ; dashed, dash-dotted and dotted curves represent contributions from partial mass ranges, as indicated 
in the legend. 



power spectrum, and A2PCF. The model based on a primordial density fluctuation field obeying xi statistics generates a 
more abundant population of massive clusters by virtue of higher probabilities for overdense regions at high redshifts. These 
collapse earlier and form cluster haloes, whereas the early quintessence model is characterized by earlier collapse owing to 
the higher linear growth factor and lower value of the critical density for collapse at high redshifts. Both models give rise to 
higher S-Z power levels, more so in the non-Gaussian model. 

We note that the two non-standard models investigated here are by no means unique; for example, lower S-Z power is 
expected in non-Gaussian models based on Xm statistics with m > 1 due to less pronounced skewness (or non-Gaussianity) of 
the density field. Likewise, early quintessence models with higher early dark energy densities are likely to further slow down 
the evolution of the linear growth factor (reducing the value of 5^) and thus give rise to higher power levels. Conversely, the 
A2PCF of S-Z clusters in the EDE model, and to a significantly higher degree in the non-Guassian model, is manifested in 
reduced correlation levels with respect to those of the standard ACDM model, in particular at high angular separations. 

EDE models with still larger quintessence densities at early times would enhance power levels and further blur the 
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distinction between such models and non-Gaussian models. On the other hand, the A2PCF provides a more adequate test 
of the viability of these models. In this regard we mention the recent work of Magliocchetti et al. (2006), who report results 
from analysis of Spitzer Space Telescope observations that presumably provide evidence for strong clustering of a galaxy 
population at z ~ 2. According to Magliocchetti et al, these observations suggest that the detected objects are very massive 
proto-spheroidal galaxies, with number densities that are considerably higher than predicted by theory. If this is indeed the 
case, it adds to a list of other observational indications of enhanced massive objects lying at relatively high rcdshifts, among 
which are the detection of structures with high velocity dispersions at redshifts ~ 4 (Miley et al. 2004) and ~ 2 (Kurk et al. 
2004), and the detection of protoclusters with masses lying in the range 2 — 9 • 10^* Mq associated with radio galaxies lying 
at z > 2 (Venemans et al. 2006). All this observational evidence points to the presence of massive objects at high rcdshifts, 
which - as shown here - may be accounted for by either a non-Gaussian, or early quintessence models. If additional evidence 
(such as reported by Magliocchetti et al.) is found for strong clustering at ^ ~ 2, it will add further support for non-Gaussian 
models. It is important to stress that these observational results pointing to strong clustering of galaxies are not necessarily 
in conflict with the lower correlations found for the non-Gaussian model presented in the previous section. This is due to 
the fact that our calculations include the contribution to the A2PCF from clusters of galaxies alone, for which, as has been 
explained above, the (non- Gaussian) bias is lower than the corresponding Gaussian bias. At lower masses, such as galaxies, the 
bias is actually stronger in non-Gaussian models, and would yield higher levels of the A2PCF. Finally, it is noteworthy that 
the linear bias factor fails to describe the actual clustering on galaxy scales, since these are already considerably non-linear. 
Models admitting non-linear corrections exist (e.g. Peacock & Dodds 1996), but are not likely to alter the general behaviour 
of the bias parameter in the non-Gaussian and EDE models. 

The non-standard cosmological models considered in this paper give rise to enchanced S-Z power by virtue of their impact 
on the cluster population. A legitimate question in this regard is whether internal cluster properties alone may be able to 
account for the apparent variance between the low as deduced from the WMAP 3-year data and the high-multipole CMB 
power excess observed by the ACBAR, CBI, and possibly BIMA experiments. We believe, however, that the answer to this 
question is negative. IC gas properties and their impact on the S-Z power spectrum have been extensively investigated in the 
literature; specifically, consequences of the gas mass fraction evolution (Majumdar 2001), isothermal, as well as polytropic 
temperature profiles (e.g. Komatsu & Seljak 2002), self-similar and entropy-driven models for the cluster core (Komatsu & 
Kitayama 1999), and the temperature-mass relation (SRS), were explored, and however significant differences were found, they 
are still incapable of settling the confiict. Recently, Roychowdhury, Ruszkowski, & Nath (2005) have explored the infiuence 
of AGN heating on groups and clusters of galaxies, motivated by the observed entropy excess in the IC gas. A consequence 
of this heating was shown to be a depletion of gas from the central regions of clusters, resulting in reduced S-Z power levels. 
Obviously, this cannot provide a remedy to the conflicting low ag and high-multipole CMB power excess. 

In conclusion we re-iterate the main result of our work: the mounting observational evidence for the presence of massive 
clusters already at relatively high rcdshifts cannot be easily reconciled with the Press & Schcchtor mass function (or variants 
thereof) in the standard ACDM cosmology with underlying Gaussian density fluctuation field. Non-Gaussian fields and early 
quintessence models may remove the inherent difficulty by allowing a more abundant population of massive clusters at high 
redshift. Moreover, such models give rise to enhanced S-Z power levels, which seem to be in better agreement with the CMB 
power excess at high multipoles observed by the CBI and ACBAR experiments. Future S-Z cluster catalogs will enable the 
construction of A2PCF from S-Z measurements, thereby providing an additional test for the viability of these models. 



5 APPENDIX 

Extension of the numerical computation of the S-Z angular power spectrum and A2PCF to EDE models necessitates a 
re-evaluation of several cosmological quantities that are affected by the presence of early dark energy. These include the 
(comoving) radial and angular diameter distances, the linear growth factor of density perturbations, the critical density for 
spherical collapse, and the overdensity at virialisation. This Appendix details the necessary modifications to the basic equations 
governing the redshift (or equivalently, the scale factor, a) dependence of these functions, and presents their numerical results. 

The EDE model studied in this work is characterised by fiat geometry, with an early quintessence density specified by 
the parameter fie, and an effective redshift-dependent coefficient of the equation of state given by 

yj(z) = — f9) 

^' 1 - [3wo log (1 + z)]/{log 1^ + log 1^) ' ^ ' 

where wo = w{z = 0). By using the redshift scaling of the density of the dark energy component corresponding to this model, 

PQ=PQoa + zf'^''^'^K (10) 
in the Friedmann equation 

(-) =— (Pm+pg), (11) 
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one obtains 

^ = H = Ho + zf + (1 - + . (12) 



1. Comoving radial and angular diameter distances 

The comoving radial distance can be obtained from equation 1121 as 

riz) = — r — (13) 

Ho Jo + z)3 +nQ{l+ 2)3[l+M-)]]l/2 ' 

with fig = 1 — flm, whereas the angular diameter distance is simply 

2. Linear growth factor 



The differential equation governing the evolution of the linear growth factor of density perturbations is 



da2 



1 flQa-^^^+'"^^^^[l + 3w{z)] + n^a' 



a 2a !r2Qa-3[i+^(^)l + n™^a-3 



dS 3 Qn 



da 2 flS nQa-3[i+M^)] + n„ 



rS^O. (15) 



3. Critical density for spherical collapse 



Using the Friedmann equations for both the spherical overdensity and the background manifold, and denoting by ata and Rta 
the scale factor and radius of the spherical overdensity at turnaround, respectively, pqt.ta, Pmb,ta, and pmc.ta the background 
density of the dark energy component, the background density of the matter component, and the collapsed density of the 
matter component at turnaround, respectively, we have 

RV SttG, , k 

-^j = —^{Pmc + PQc) - ^ (16) 
dV SttG 



g (pmb+PQb), (17) 

where the curvature k is assumed to vanish in our model, and the indexes c and b refer to collapsed and background components, 
respectively, the following differential equation may be obtained: 

.dx) a:-i + i/a;-iat^'""'V(a;-ata)3»(--*») 

with X = a/ata, y = R/Rta, v = pQb,ta/ pmb,ta, and ( = Pnic,ta/pmb,ta = {Rta / CLtaY'^ ■ Siuce X = and y = Q sX t = Q, whereas 
X = 1 and J/ = 1 at turnaround, it is possible to recast equation [18] in the form 

dy 

1/2 

(19) 



« [C(y-' - 1) + v{y-^{ataC-^/'3Y-^(a,^C-'/^)/(^y . („,„^-l/3))3-(y.(a,„C-l/=')) _ 1)] 

dx 



\x'^ + i/ar-ia?™'"'"' /{x ■ ata)3-{--*")] 



1/2 ■ 



where we have used the relation Rta = ataC ■ Equation (|19p can be solved numerically to yield the sought value of C,, and 
the critical density for spherical collapse can then be found as (Zeng & Gao 2005a): 

5. = \at^ [C' + i^ataC^'') D{a,), (20) 
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where D{ac) is the hnear growth factor at the time of coUapse. 



4. Overdensity at virialisation 



The overdensity of the viriahsed matter component can be found by relating the background and coUapsed matter 
densities at turnaround with their counterparts at the time of collapse: 

^ _ Pme,c _ pmc,ta {Rta / Rc)^ _ ^ f \ ^ ^^^^ 

pmb,c pmb.ta {flta / O-c)'^ ' \CLta J 

or if we are rather interested in the virialised overdensity in terms of the critical density at the time of collapse 

a; = = Qmia.) ■ A„, (22) 

pc, c 

where we have defined A = Rc/Rta- It remains to specify this parameter. This can be done by requiring energy conservation at 
turnaround and virialisation, in conjunction with the virial theorem applied to the virialised phase of the collapsed overdensity. 
Before proceeding along this line, we remark that there is an ongoing debate regarding the assumption of energy conservation 
within the collapsing halo in models with dark energy featuring w ^ —1 (e.g. Basilakos & Voglis 2007, Zeng & Gao 2005b, 
Maor & Lahav 2005). In fact, Maor & Lahav show that in such models the assumption of energy conservation breaks down. 
Here we choose to ignore this issue, and use energy conservation as a zero order approximation. It is also of interest to note 
that according to Wang (2006), under certain conditions (namely, for u/d^ < 0.01), the problematics associated with energy 
conservation do not affect the virialisation process. Although this factor is somewhat larger in our adopted model, we believe 
that our results are not significantly affected by this ambiguity. 

To proceed, the potential energy stored in the matter component of the spherical overdensity is the usual 

Ug = (23) 

whereas the potential energy associated with the dark energy in the spherical overdensity is (Horellou & Berge 2005) 

[/g = (l + 3^)pg.^i?^ (24) 

At turnaround the spherical overdensity has only potential energy, while at virialisation it has both kinetic and potential 
energies: 

{T + U)ta = {T+U),. (25) 
Using the virial theorem, Tc — —1/2Ug + Uq, so 

^Ug,c + 2Uq,c = UG,ta + UQ,ta. (26) 

Substitution of equations (|23|) and (|24p in equation (I26p (and performing some algebraic manipulations) leads to the following 
equation: 

Rta ^ 47r[l + m{ata)]pQMR^a " 6M 

Rc 87r[l + 3«J(ac)]pQ,ci?3 - 3M ' ' 
Using the fact that the mass can be represented by either M = 'iTT/3pmc,taRta or M — Aiv/Spmc.cRc, this can be written as 

R.a + 

Rc 



2[l + 3-(ac)]i^(ff) {^) i 



(28) 



- 1 



Given the collapse scale factor Oc, or, alternatively, the collapse redshift Zc, all of the density parameters appearing in this 
equation can be found. Substituting A = RcRta, the last equation can be written as 

m - 2 1 



2n2A3-l A' 



(29) 



2n2A^ - (nl - 2)A- 1 = 0, (30) 

where ni = [1 + 3w{ata)]C~^ i^Q,ta/^rnb.ta), and n2 = [1 + 3in{ac)](~^ {ilQ,cata) / i^mb,cal) can aU be evaluated given the 
collapse redshift. The solution for A should obviously lie in the range ^ A ^ 1. 
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Figure 5. Numerical results for the linear growth factor, critical density for spherical collapse, and the overdensity at virialisation are 
illustrated in the upper left-hand, upper right-hand, and bottom panels, respectively, for models I (continuous) and II (dashed). 

Numerical results for the linear growth factor, critical overdensity for collapse, and overdensity at virialisation in the 
EDE model are compared with the corresponding results for the standard ACDM model in Fig. [S] 
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